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In the first moments of a relativistic heavy ion collision explosive collective flow begins to grow 
before the matter has yet equilibrated. Here it is found that as long as the stress-energy tensor 
is traceless, early flow is independent of whether the matter is composed of fields or particles, 
equilibrated or not, or whether the stress-energy tensor is isotropic. This eliminates much of the 
uncertainty in modeling early stages of a collision. 
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O . I. INTRODUCTION AND BASIC THEORY 

o 

At ultra-relativistic energies, heavy ions collide and produce enormous energy densities, exceeding several 
(j ' GeV/fm 3 , well above the threshold for dissolving hadrons into partons and melting the QCD vacuum condensates 
• Due to partial transparency the deposited energy is borne with a large collective velocity gradient along the 
beam axis (which will be labeled the z axis). In the limit of infinite beam energy, the collective expansion along 
the z axis is boost invariant [l|, with v z = z/t, which corresponds to zero acceleration in the z direction. Since 
the matter has no initial collective velocity in the transverse direction, transverse expansion is pressure driven 
with the transverse collective velocities growing until break-up at which time particles free-stream toward the 
detector. Several clear signals of the collective nature of the transverse expansion have been observed. First, the 
average transverse energy of protons, kaons and pions are ordered by mass [3[, as expected since the collective 
velocity adds more energy to a more massive particle. Secondly, for non-central collisions the initial transverse 
energy density profile is anisotropic, which leads to anisotropic transverse pressure gradients and thus anisotropic 
£h ■ flow The observable, 

v 2 = i(cos2(>-&.p.)), (1) 

>: 

, quantifies the anisotropy, where 4> T .p. is the angle of the reaction plane. Observed values of v 2 reach into the tens 
O^l i of percent [J] and are consistent with expectations of hydrodynamic flow [f|. Finally, two particle correlations, 
which are measured as functions of the two momenta p\ and p 2 , provide six-dimensional femtoscopic pictures of 
outgoing phase space clouds. The six-dimensional structure has numerous features expected from boost-invariant 
flow along the z axis and from strong transverse flow Q. 

Comparing data to simple models reveals solid evidence for the existence of strong collective flow. However, 
detailed comparison with full dynamic models are required to infer quantitative information about the equation 
of state, viscosity, or other properties of QCD matter. Modeling relativistic heavy ion collisions is complicated 
since the matter traverses three distinct stages during the ~ 15 fm/c expansion, each of which involves different 
degrees of freedom. In the final 5-10 fm/c of the collision, thermal equilibrium is lost and hydrodynamics, 
even viscous hydrodynamics, is inapplicable. Collisions at this point are binary and microscopic simulations 

In the intermediate stage, local equilibrium 



OO 

o 



5-H 



or Boltzmann codes can be applied with some confidence [7 



is sufficiently maintained to justify viscous hydrodynamics [9|, JJ. This is fortunate, since the quantities of 
greatest interest - temperature, pressure, and viscosity - form the basis of the description. It is during this stage 
that the majority of collective flow develops. The first stage, times less than < 1 fm/c, is the most uncertain and 
theoretically contentious. Descriptions might be based on classical QCD fields [ll|, 0> or on partons [HI]. The 
partons might be highly non-equilibrated with a significant fraction of the energy tied up in high-energy jets [Til ]. 
The initial transverse acceleration is driven by the conservation laws of the stress-energy tensor, d a T a ° = 0, and 
given the very different stress-energy tensors implied by the proposed models, one might expect very different 
flows to be generated during the early stage. The importance of understanding the early phase is underscored 
by noting that although the pre- hydrodynamic stage lasts only ~ 1 fm/c, the acceleration during that stage has 
a relatively higher impact on the evolution of the collision. This is for the same reason that the start plays a 
relatively large role in a 100-m sprint. 

The purpose of this article is to explain a non-intuitive result - that the impulse to the collective flow provided 
by the initial stage is identical for a large class of simple models, even though the pictures yield very different 
evolutions of the stress energy tensor. The requirements for the universal behavior are threefold: 

1. Longitudinal flow has a boost-invariant profile, v z = z/t. 

2. The stress-energy tensor is traceless. 
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Model 




T Z2 


Too 


Longitudinal electric field 


Too 


—Too 


~ const ant 


Free streaming massless particles, 
two-dimensional relativistic gas 
or fields from incoherent longitudinal currents 


Too/2 





~ l/t 


Ideal hydrodynamics of massless gas 


Too/3 


Too/3 


~ iA 4/3 



TABLE I: Elements of the stress energy tensor for simple models at early times. A wide variety of reasonable values for 
the stiffness of the transverse equation of state are represented by the models. 



3. The anisotropy of the spatial components of the stress energy tensor is independent of the transverse 
coordinate and depends only on the Bjorken time r = \Jt 2 — z 2 . 

The first assumption is reasonable for 100A GeV collisions at the Relativistic Heavy Ion Collider (RHIC) 
and very well justified for the very high energy heavy ion collisions to be performed at the LHC. The second 
assumption is well satisfied by any description based on massless particles or weakly interacting gauge fields. 
The final assumption is more subtle: for instance, if a system suddenly changes from longitudinal fields to 
thcrmalizcd particles, then the anisotropy of the stress-energy tensor will also change suddenly. As a measure 
of the anistropy, one can define the quantity, 

k = (T xx + T yy )/2T 00 , (2) 

and if k varies mainly with r, rather than with the transverse coordinates, the final criterion for universality is 
satisfied. Once the criteria for universal behavior are satisfied, the hydrodynamic stage can be initiated with 
known flow fields, which should depend only on the shape of the initial energy density profile. The simple 
expression for the flow fields derived here reduces, if not eliminates, the need for detailed modeling of the initial 
stage, and isolates the final evolution from contentious issues concerning the pre-equilibrium stage. 
Collective flow is driven by the stress energy tensor, 
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Txo T xx 
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V T Z2 



(3) 



where it has been assumed that one is modeling only along the y — z — axis, so that v z = (although 
d z v z = l/t), and x refers to the radial direction with v y = by reflection symmetry about the yz plane. 
For central collisions, the initial stress-energy tensor has T xx = T yy and is diagonal, with T xx , T yy and T zz 
representing the force per unit area pushing in the three directions. 

If one assumes a small velocity v x and neglects the longitudinal expansion, one finds a simple expression for 
the radial acceleration, 

dv x = dT xx /dx ^ 
dt T 00 + T XX 

This immediately gives the impression that the transverse acceleration is driven by T xx , which differs substan- 
tially for different pictures of the initial stage listed in Table [IJ However, the inclusion of longitudinal expansion 
significantly alters the expression for the radial acceleration. To see this we consider the full expression for the 
conservation of momentum current, 

dtTo x — —d x T xx — d z T xz — d y T xy , (5) 

Expressions for T xz and T xy for small z and y can be expressed using boost and rotational invariance, 

z y 
T xz {x, y = 0, z) ps -Tox, T xy (x,y, z = 0) ps ~{T X x - T yy ). (6) 



Inserting these into the last two terms in Eq. ([5]), 



Toa: {T xx Tyy ) 



d t T 0x = -d x T xx - - ^ ^ . (7) 

t x 

Similarly, conservation of the energy current, 

dt Too = —d x To x — d y To y — 8 z Tq z , (8) 
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can be rewritten using boost and rotational invariance to re-express Tq v and Tqz, 



Toz(x,V = 0,2) » (Too + T zz )j, T Oy (x,y,0) ps ~Tq x , (9) 



which transforms Eq. (jSJ) along the y = z = axis to 



Combining Eq.s and JTDJ, 



5tT 00 = - i (Too + T zz ) - ^ + i ) 7i ,, . ( 1 1 > ! 

00/ -too -too 1 



d/ -~ ux \ ~x-^xx , ^-zz^-vx 

-tOO / -tOO -innt 



Only early times are considered so terms of higher order in Tq x /T 00 and the term proportional to (T xx — T yy ) 
are discarded. Since Tq x /Tqo rises linearly in time, one can solve for the effective acceleration a{x) defined by, 

T 

— - = a(x)t, for small t, (12) 



Too 
a(x) 



r) T 

^x-t X 



Too — T zz 



If the stress energy tensor is traceless and if T xx = T yy = «(i)Too, that is the anisotropy of Ty does not depend 
on x, one obtains, 

«(*) = ^ 

-t XX T" -t yy 

d x Too 
2Tqo 

which is independent of k. The principal result of this study, 

Tq x ^ dxToo 
Too 2Tnn 



t, (14) 



is surprising in that the acceleration is independent of the anisotropy of the equation of state. The extra strength 
of the transverse acceleration due to longitudinal flow is also remarkable. For the case of ideal hydrodynamics, 
Tqx/Too increases 50% faster than what one would expect from Eq. ((4]), which would have been applied if there 
were no longitudinal expansion. 

The ratio Tq x /Tqq can determine velocity once one has made an assumption about the equation of state and 
viscosity to be applied in the hydrodynamic stage. Combining the constraint of Eq. (fl"4"|) with the Navier-Stokes 
equation and knowing the initial profile for Too is sufficient to determine all of T"' 3 '. However, the success of this 
procedure depends on whether the universal flow persists to sufficiently late times so that the velocity gradient, 
dv z jdz ~ has subsided to the point that the Navier-Stokes conditions for the stress-energy tensor can be 
applied. Estimates for such a condition are close to < 1 fm/c. The goal of the next section is to exactly calculate 
Tox/Tqo for several simple models, with widely different anisotropies of Ty, to see whether the flow remains 
universal out to times near 1 fm/c. 

Although Tq x /Toq has a universal behavior for the models being considered here, other measures of flow will 
vary between models. For example, one might use the collective velocity, u, of the matter, i.e., the velocity 
for which an observer would measure Tqi = 0. A common measure of elliptic flow [lol [l5| uses the spatial 
components of the stress-energy tensor, 

£ = / dxd V (T™ ~ Tyy) ^ 
p J dxdy (T xx +T yy Y 

which for free particles would be proportional to v%. All these measures are defined through the stress-energy 
tensor. The quantity representing the development of collective flow should preferable be chosen such that it 
is preserved during a rapid change of the microscopic degrees of freedom. Examples of such a change might be 
the sudden decay of longitudinal fields or the rapid isotropization of matter. During such a transition, the basic 
conservation of the stress-energy tensor should remain valid, even as any of its ten independent components 
are allowed to change suddenly. One can consider a hyper-surface for which the microscopic form of matter 



4 



is discontinuous on one side to the other. For example, for one side of the hyper-surface one might assume 
longitudinal fields, while applying ideal hydrodynamics on the opposite side. At any point on the hypersurface, 
one can find the four- vector n a which is orthogonal to the surface. If n 2 = +1 it is possible to boost to the 
frame where n = (1, 0, 0, 0). In this frame, neighboring points all undergo the transition simultaneously and the 
only discontinuity is in the time direction. If one integrates the conservation equations across an infinitesimal 
time element, one finds the four constraints 

rt+St 

= / (d t T 0a + diT ia ) (16) 

Jt-St 

= T aQ (x, y,z,t + St) - T a0 (x, y,z,t- 5t). 

If the hyper-surface were defined by a unit four- vector where n 2 — —1, one could boost to a frame where the 
discontinuity was locally static. In that case one would integrate the conservation equations across the surface, 
and if the discontinuity were in the z direction, one would find: 

= T az {x, y,z + Sz, t) - T az {x, y,z- 5z, t). (17) 

These are the Rankine-Hugoniot conditions for a shock wave. For either case, energy-momentum conservation 
across a discontinuity is stated as 

n a T a(3 {x+) = n a T a/3 {x-). (18) 

The sudden changes of state most often imposed in relativistic heavy-ion collisions, e.g., changing from fields 
to particles, tend to be invoked at a constant Bjorken time, r = \/t 2 — z 2 , which corresponds to n — (1, 0, 0, 0) 
for z = 0. This would force T 00 and T oi to remain constant across the discontinuity, and the flow as defined by 
Toi/Too would be unchanged. In contrast, since both the flow velocity u and the elliptic anisotropy e p require 
using the spatial components of the stress-energy tensor for their defnition, they could change instantaneously 
without violating energy-momentum conservation. Thus, these measures can provide misleading insight into 
whether flow has developed. For instance, if one considers free streaming particles, e p remains zero. But if those 
particles instantaneously thermalize at a constant Bjorken r, the stress-energy tensor will instantaneously adjust 
itself while preserving the four spatial components To a . At this point, e p suddenly becomes finite. Furthermore, 
to the degree that Toi /Too is universal, the collective flow will be determined by the new state of matter rather 
than the previous state. For example, if one believes the matter becomes suddenly hydrodynamic at r = 1 fm/c, 
all the models listed in Table U would lead to the same flow for r > 1 fm/c, with the same collective velocities 
u and the same e p . 

In the next section, simple examples are presented where the universality of flow is tested for sudden transitions 
to hydrodynamics. Similar universal behavior would apply if the transition were to viscous hydrodynamics, 
though the collective velocities u and the elliptic anisotropy e p would be different than for transitions to ideal 
hydrodynamics. Thus, once one has chosen a time To at which to begin the hydrodynamic prescription, along 
with the viscosity coefficient, the collective velocity at To is then independent of the properties of the matter 
assumed for t < to, as long as the three criteria listed earlier are satisfied. 

II. THREE MODELS 

In the previous section it was shown that for small times one finds universal behavior of the flow as defined 
by Tq x /Too- In order to assess whether universal behavior applies for times of the order r < 1.0 fm/c, exact 
solutions are found for three different models: 

• Ideal hydrodynamics of an ultra-relativistic gas with P = e/3. 

• Coherent non-interacting electro-magnetic fields (field points in the same direction initially) that follow 
Maxwell's equations. This can be thought of as the field generated by two oppositely charged capacitor 
plates receding at ±c. 

• Incoherent electromagnetic fields generated by a random ensemble of ± charges receding from one another 
at ±c. 



Each of these cases involves a traceless stress-energy tensor, 

TxX H~ Tyy -\- T zz 6, 



(19) 



5 



but with very different asymmetries between the initial transverse pressure T xx — T yy and the longitudinal 
pressure T zz . The measure of this initial anisotropy, k, which is the ratio of the transverse pressure T xx to the 
energy density e as defined in Eq. ((2]), is 1/3 for the hydrodynamic model, 1/2 for the incoherent fields and 
unity for the coherent fields. Non-interacting massless particles would also result in n — 1/2. Colliding particles 
with insufficient cross section to thermalize would have 1/3 < k < 1/2. Thus, these three models seem to well 
span the range of possibilities. For the solutions presented in the next section, each model is initialized with 
the same transverse energy density profile and with zero collective velocity. 

The first model is ideal hydrodynamics with a simple equation of state P = e/3. For ideal hydrodynamics 
one assumes that the stress-energy tensor has the form, 

T afi = (P + e)u a u p - Pg ap . (20) 

For this model T xx — T yy — Too/3 at early times, i.e., k = 1/3. The model is solved numerically, beginning at 
t = 0. Solving the equation of motion, d a T a @ = 0, is complicated by the fact that the energy density is singular 
at r = 0. For that reason, the equations of motion are manipulated so that one solves for the quantities 

» = ijH (2i) 

J 00 

U = Toor 4 / 3 . 

The equations of motion for w and U are straight-forward to generate from Eq.s and (|10|). with the advantage 
being that the quantities have no singular behavior. 

For the second model we consider a coherent electromagnetic field generated from two oppositely charged 
capacitor plates receding from one another at ±c. This is not particularly physical, as the fields are coherently 
pointing along the same direction, whereas in a heavy-ion collision the chromo-electric fields would be random 
with a coherence length set by the saturation scale flU [l7j . This model is chosen because it has the extreme 
T xx ~ Too as the initial equation of state. For a single pair of opposite charges originating from x = y = 0, 
Lienart Wiechart potentials can be used to generate the electric and magnetic fields, 

/>OC 

A z (r,t) = 2q / dt'S(x 2 + y 2 - (t - t') 2 ), (22) 
Jo 

E s (r,t) = 4qS(r 2 -t 2 ), 
B$(r,t) = -E z {r,t). 

For receeding plates with charge densities p(x,y), one can integrate over the charges and find, 

E z (x,y,t) — 2 J d<f> p(x ~ t cos <fi, y — t sin 0), (23) 
B x (x,y,t) = 2j d(f> p(x — t cos <j>, y — t sin 0) sin (j>, 
B y (x,y,t) = —2 J d(f) p(x — t cos 4>, y — t sin</>) cos cj). 

The fields can then be used to generate the stress-energy tensor, 

(24) 



Tqq — 




B 2 X + 


Bl), 


T — 




B 2 X + 


B% 


T — 
± yy 




B 2 + 


B% 


T — 


\(-E\ 


+ Bl 




Tq x = 


—E z B y , 






T y = 


E Z B X - 







The initial transverse distribution of the electric field matches the profile of the charge density, and since 
Too ~ E 2 Z initially, the gaussian radius characterizing the electric field is larger than that of the resulting energy 
density profile by a factor of \[2. The charge density was then chosen with the larger radius so that the resulting 
energy density profile would be the same as the other two models. 



G 



The third model is also based on the evolution of classical fields, but assumes that the fields resulted from 
a distribution of point charges ±g, receeding at ±c. The resulting fields from a single point particle are easily 
calculated from Eq. (|22|) . Although the fields have a random sign, the resulting stress-energy tensor is always 
positive and the energy density always moves outward. Thus, one calculates the stress-energy tensors for a 
single point-charge pair, then integrates over the density of such charges to find the stress-energy tensor. For 
such a case the transverse magnetic field and the longitudinal electric fields have the same strength, which 
makes T zz — and k — 1/2. Effectively, each point sends out an electromagnetic pulse which behaves exactly 
the same as massless partons being emitted from a point source at z = t = 0. ft is not surprising that the value 
of k = 1/2 is identical to what one would obtain from a non-interacting parton picture, where all the partons 
were emitted at z = t = 0. 

Complete incoherence, as assumed in the third model, is unphysical as it assumes perfectly point-like charges. 
For a beam energy E, the uncertainty principle precludes assigning the starting points for z and t to a region 
less than ~ 1/E. The transverse sizes of the individual charges is also limited by whatever radiative dynamics 
are responsible for the creation of the charge exchanges. In QCD this coherence length is referred to as the 
saturation scale [HI]. For finite-sized charges of correlation length A, the stress-energy tensor behaves like the 
coherent limit for r << A and like the incoherent limit for r >> A. From the delta function forms for the electric 
and magnetic field in Eq. (|22[) one can see that the energy density, which involves squaring the fields which 
are described by delta functions, is infinite in the limit of A — * 0. In fact the radiated density per unit rapidity 
behaves ~ a/ A. While it is clear that a realistic evolution of the fields is more complicated than either the 
coherent or the incoherent limit even in the limit of little interaction, if one considers an ensemble of randomly 
charged small but finite-sized Gaussian packets, k would still depend mainly on r, and not x or y. The three 
criteria for universal early flow would still be satisfied. 



III. RESULTS 



Here, we present the evolution of Tq x /Toq for the three simple models described in Table |TJ each with widely 
different initial anisotropies for T%j. For each model, we assume a Gaussian energy density profile, 

T 00 (z, y )c<expj-^-2lf}, (25) 

with radii R x = R y = 3 fm for the calculations used to investigate radial flow. The transverse velocity and flow, 
u x and To x /Too respectively, are calculated as a function of the transverse coordinate x along the y = z = 
axis for three early times: 0.3, 0.6 and 1.0 fm/c. Additionally, one can calculate the collective velocity that 
would ensue if the system were to suddenly change at those three times into pure hydrodynamic flow with a 
simple equation of state P = e/3. These are found by comparing to the form for the stress energy tensor for 
ideal hydrodynamics as defined by Eq. (j2"D|) . One solves for u' x and e by matching to Too and Tq x , 

4 2 1 4 
Too = 2 eu 'o ~ 3 e ' = 3 eu a u 'x- (26) 

The velocity u' x will depend only on the ratio Tq x /Toq. Thus, if Tq x /Tqo exhibits universal behavior, so will u' x . 
For the hydrodynamic model v! = u, whereas they differ for the other two models. From Eq. (fi~4")) . one can 
exactly calculate To x /Too to lowest order in the proper time r, 

Too 2R X 

Using Eq. (|26[) one can also find u x and u' x to lowest order in r, 

XT 

- - (iT^r < 28 > 

, XT 



(4/3)ig" 



Figure [T] displays u x , T 0x /T 00 and u' x for all three models as a function of x for three times, 0.3, 0.6 and 1.0 
fm/c. As expected, Tq x /Tqq and u' x exhibit nearly universal behavior, with differences of a few percent by 1 fm/c. 
The linear-in-time approximation, given in Eq. (|27p . is also displayed for Tq x /Tqq. The linear approximation is 
remarkably effective for the first 1.0 fm/c. As expected from Eq. (|2"5|) . the flow velocities, u x , differ. What is 
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FIG. 1: (color online) Lower Panel: The collective velocity profile is displayed for three models at three different times, 
0.3, 0.6 and 1.0 fm/c. Ideal hydrodynamics (green triangles) has the greatest transverse radial collective flow even though 
it had the smallest transverse pressure, T xx , of all three models. The evolution of a non-interacting coherent electric field 
(blue squares) had the highest pressure, but the smallest flow. Calculations based on electromagnetic fields arising from 
incoherent currents (red circles) would be the same as for non-interacting partons. 

Middle Panel: The flow ratio Tox/Too is nearly universal for all three models. The symbols are the same as in the lower 
panel. The solid lines represent the linear approximation, ~ r, given in Eq. (|28[) . 

Upper Panel: The collective velocity assuming that the matter suddenly behaves as if it were ideal hydrodynamics at 
the prescribed time. Since this ratio depends on Tq x /Tqq, it is also nearly universal. 



surprising is that the models with higher values of k lead to lower velocities. This is opposite to the trend one 
would obtain if there were no longitudinal flow and Eq. Q would be have been applicable. 

Elliptic flow was evaluated by considering emission from an initial energy profile characterized by R x = 2, 
and R y = 3. As a measure of elliptic flow, e p , defined in Eq. (fT5]) . is calculated for the two models based 
on coherent and incoherent fields. Results for the hydrodynamic model are skipped because that model was 
predicated on radial symmetry, althou gh c alculations have been done previously for ideal hydrodynamics with 
and without transverse thermalization [19| . Assuming a sudden transformation to ideal hydrodynamics at r, e' p 
was also calculated using the same method to calculate u' x used for the radial case above. Figure O shows both 
e p and e' p as a function of r for both models and compares them to the small-r expansion, e p ~ r 2 . The small-r 
limit is found by calculating u x and u y for small times from Eq. (|28[l for the hydrodynamic model, k — 1/3. 
The collective velocities are then 



u 



(hydro) M 3£^ u ( hydro ) w 3^ (2Q) 



Using Eq. (|20|) for the stress-energy tensor, one can then calculate the elliptic anisotropy with some straight- 
forward integrals of Gaussians, 



8 




0.0 0.5 1.0 

r (fm/c) 

FIG. 2: (color online) Lower Panel: The elliptic anisotropy e p as defined by Eq. (|15[) for the case of coherent non- 
interacting initially-longitudinal electric fields (blue squares) and for fields driven by incoherent currents (red circles). 
The incoherent case yields e p is zero, exactly as one would obtain with non-interacting particles. The solid line shows 
the lowest-order (in r) quadratic contribution for ideal hydrodynamics. 

Upper Panel: Assuming that the matter suddenly behaves as if it were ideal hydrodynamics at time r, e' p represents the 
anisotropy of the altered stress-energy tensor. The result is close to the quadratic form approximating the behavior of 
ideal hydrodynamics. 



As expected, the two models agree with this simple quadratic form for e' p , but differ very substantially for e p . 
In fact, for the model with incoherent fields, e p remains zero for all times. This follows from the fact that 
each point source contributes incoherently to the stress-energy tensor, and each point source has zero elliptic 
anisotropy. 

Once a system has decoupled, the anisotropy e p can be equated with the angular anisotropy V2, 

(v 2 ) ~ ^e P , (31) 

where the average (• • • ) refer to an average over particles in a central rapidity bin weighted by pl/m t . Even 
though the values of e p are < 10% of the v 2 observed experimentally, the contribution from the first fm/c 
is substantial. Since e p grows quadratically in time, it is important to generate a rate of change, de p /dT, as 
quickly as possible. The first fm/c is especially important in elliptic flow analyses for two reasons. First, one is 
considering non-central collisions which are smaller in overall volume and thus of shorter duration, and second, 
elliptic flow saturates earlier than radial flow [l^| . 



IV. SUMMARY 



The existence of universal flow patterns for the first < 1.0 fm/c of a relativistic heavy ion collision has a pro- 
found impact on the modeling and interpretation of heavy ion collisions. It eliminates many of the uncertainties 
plaguing the pre-thermalized stage. For example, if one were to use viscous hydrodynamics beginning at r = 1 
fm/c, the initial profile for Tqi/Tqq would be determined by the universal conditions shown here. Given that 
the contribution to the final-state flow from the first fm/c could be of the order of 10-20%, it makes detailed 
modeling of the pre-thermalized stage unnecessary if one is only interested in the development of the evolution 
of the stress-energy tensor at later times. 

This does not, by any means, make theoretical investigations of the pre-thermalized stage irrelevant. Uncer- 
tainties in the shape of the initial profile would remain, including questions about the magnitude of the initial 
energy density and the microscopic structure. Even though two pictures result in the same flow fields, they 
might have very different microscopic structure. Differing densities of quarks, gluons and kinetic temperatures 
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should affect a variety of other observables such as electromagnetic probes [20|, |2lj, l22|, jet quenching [23[, or 
charge balance functions 0, [H| . 
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